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STATISTICAL THEORY OF DROUGHTS 


E. J. Gumbel 
Adjunct Professor, Columbia University, New York* 


ABSTRACT 


Droughts are analyzed by the asymptotic theory of smallest values of a 
limited statistical variate. The extremal probability paper used for floods 
is used for the logarithms of droughts. If the lower limit of the discharges 
is assumed to be zero the probability function of the droughts becomes a 
straight line. If the lower limit exceeds zero, the three parameters in the 
probability function are estimated by the method of moments. A given statis- 
tical criterion indicates whether or not the lower limit may be assumed to be 
zero. Observations on the droughts of 13 rivers analyzed by this procedure 
show a close conformance with the theory. Therefore, extrapolation is per- 
missible if the basic conditions prevail. 


1. Introduction 


The asymptotic theory of extreme values taken from an unlimited initial 
distribution of the exponential type is now widely used in the analysis of floods 


[2, 3, 9, 10].** The fact that this theory as well as many others assumes no 
upper limit, although shocking to some engineers, has not led to major diffi- 
culties. The main difficulty encountered has been the occurrence in smaller 
rivers of extraordinary floods exceeding the magnitude expected for the given 
number of observations. 

The theory of extreme values is also an appropriate tool for the analysis 
of droughts which are defined as the annual minima of discharges. Up to now, 
the theory of smallest values taken from an unlimited distribution of the ex- 
ponential type has been used [12]. The droughts are plotted in decreasing 
order on probability paper designed for the theory of extreme values. How- 
ever, this procedure leads to a sort of discontinuity. Instead of one line, two 
lines have to be used, one dealing with moderate, and the other with severe 
droughts. The first part is interpreted as if it did not belong to the extreme 
values proper, but still to the initial distribution. Only the latter part serves 
for extrapolation. This leads to a loss of about 37% of the information fur- 
nished by the observation. 

In the following another asymptotic theory of smallest values is used. It 
takes into account the fundamental difference between floods and droughts. 
For floods the lower limit of the discharges is of no practical interest. It does 
not affect the theory because the same asymptotic probability of largest values 
holds, as long as the initial variate is positive or unlimited to the left. How- 
ever, for droughts being smallest values, the lower limit is assumed to be 


*Work done in part as Consultant to Stanford University under the sponsor- 
ship of the Office of Naval Research and in part under grant from the Higgins 


foundation. 
**The numbers in brackets [ ] refer to the bibliography. 
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(1) zero, or (2) a small positive value €. Under adequate climatologic condi- 
tions, the first assumption may hold for small streams. The second assump- 
tion must be ‘ntroduced for large rivers. 

To use the lower limit, consider the three asymptotic probabilities of ex- 
treme values given in Table I where ®(x) (and P(x)) is the probability that a 
largest (smallest) value is equal to or smaller (larger) than x. The second 
asymptotic probability was constructed by Fréchet [4]. The other two types 
ar: due to Fisher and Tippett [5]. A systematic theory linking the three 
asymptotes to the corresponding types of initial distributions was given by 
R. von Mises [7]. 

Probability tables for the first asymptotic distribution of extreme values 
have recently been published by the National Bureau of Standards [8]. The 
first and the third asymptotic distributions of extreme values are related by 
a logarithmic transformation which will be used in the following. 

All three types involve a location parameter u and a scale parameter 1/a. 
In addition the third (limited) type involves an upper limit or a lower limit 
€. For x = u, all six probabilities are 


(1.1) ®(u) = 1/e = P(u) 


Up to now, the first type has been used exclusively. In the following, the third 
type is used. Its analytic properties are investigated and its averages and 
moments calculated. Then the first asymptotic probability of largest values 
(used for floods) and the third asymptotic probability of smallest values (used 
for droughts) are compared and their similarities noted. Two methods for 
estimating the parameters are given corresponding to the cases (1) « = 0 and 
(2)€ >0. Finally, the theory is compared to observed droughts. 


2. Lower Limit Zero 


Let x be a drought, measured in cubic feet per second (cfs) or cubic feet 
per second per square mile, or cubic meters per second. Then P 3)(x) is 
the probability of a drought exceeding x. Suppose the lower limit is zero, then 
this probability is, from Table I, dropping the index (3) 


(2.1) P(x) = exp [-(x/u)@] 


where P(O) = 1 means that all droughts exceed the discharge zero. The value 
x = u, the drought which will be exceeded 36.788 %of the time may serve to 
characterize a given river. Therefore it is called the characteristic drought. 
The median drought x is 


(2.2) = u (Ig 2)1/a 
where lg stands for the natural logarithm. The modal drought X obtained 
after two differentiations of (2.1). 
(2.3) ® = u(1 - 1/a)1/4 
is smaller than the characteristic drought u. A mode exists only for 1/a< 1, 
and the mode 

precedes 


> 
(2.4) equals the median if 1/a ( = 0.30685 


exceeds < 


Hence, the existence of a mode and its location relative to the median depend 
on the scale parameter 1/a. In the case 1/a = 0.30685 the distribution is 
nearly symmetrical but not normal. 
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For numerical use of (2.1), the transformation 


(2.5) x=e7; u=eY 


is introduced. Then the Frobability /7 (z) that the new variable exceeds z 
becomes 
(2.6) IT(z) = exp [-ew@-v)] 


Table I shows that this expression is the first asymptotic probability of 
smallest values and depends only upon the reduced variate 


(2.7) -y = a@' (log x - log u) 
Here, log stands for the common logarithm, and @ and @' are related by 
(2.7') a' = 2.302594; 1/a' = 0.43429/a 


The reduced variate y is a linear function of log x. Therefore the common 
logarithms of droughts can be treated by the methods used for floods. In 
particular, extremal probability paper can also be used for droughts if their 
common logarithms are plotted on the linear (vertical) scale. The common 
logarithm of the mth drought (in decreasing order) is plotted at the frequency 
m/(N+1) where N stands for the total number of observed years. If the same 
drought x is observed k times, say, a mean rank 


(2.8) m = Vm(m +k - 1) 


is assigned to x and is then plotted at m/(N+1). This procedure is used for 
all droughts except the smallest. If the smallest is observed several times, 
all values are plotted since the number of observations should be preserved 
on the graph. Graphs 1 and 2 show droughts of the Colorado and several other 
rivers traced on logarithmic extremal probability paper, prepared by the 
University of Michigan’s School of Public Health. To facilitate reading, the 
probabilities P(x) and return periods T(x) are traced on horizontal parallel 
scales. The return period T(x) of a drought more severe than x is defined as 
usually by 


T (x) = 1/[1 - P 
whence from (2.1) 
(2.9) T (x) = 1/(1 - exp [-(x/u)®]) 


The droughts decrease in amount and increase in severity with increasing re- 
turn periods. 

The procedure for droughts is thus strictly analogous to the procedure for 
floods. The interpretation and the influencé of the parameters, however, differ 
in the two cases. For floods, the mode is u and exists independently of the 
value of the other parameter. The mode precedes the median and the skew- 
ness is fixed. For droughts the mode and the median given by (2.3) and (2.2) 
and the skewness depend upon 1/a. For floods this parameter has the dimen- 
sion of a discharge; for droughts it is dimensionless. Therefore the slope of 
the line (2.7) is independent of the measures chosen for the droughts. 

Consider now the asymptotic behavior of severe droughts. As shown pre- 
viously [6] the natural logarithm of the return period traced on extremal proba- 
bility paper converges toward the reduced variate y. Therefore the relation 
(2.9) can be written asymptotically for severe droughts 


(2.9') - lgT ~ a(lgx - lgu) 
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and log x converges to a linear function of log T while, for floods, x itself 
converges to a linear function of log T. The parameter 1/a is the negative 
of the derivative of the logarithm of the drought with respect to the logarithm 
of the return period. Therefore drougnts decrease asymptotically as a power 
of the return period, while floods increase as the logarithm of the return 
period. The return period of the zero drought is infinity. Consequently, the 
assumption of a lower drought limit of zero does not mean that a river may 
be expected to be dry at a certain date. 

The drought as a function of the return period T may be written asymp- 
totically, from (2.9') for severe droughts, as 


(2.10) 


Therefore the drought Xor corresponding to the return period 2T is related 
to the drought x, by 


(2.11) xoT ~xq( 


If return period is doubled, the expected severest drought is asymptotically 
the original drought divided by 21/a while the doubling of the return period 
leads only to the addition of (lg2)/a to a flood. 

Conversely, consider the return period T' during whichadrought decreases 
to one half its initial amount. From 


& 
and (2.10) it follows that this return period T' is asymptotically 2“ T where 
T is the original return period. On the other hand the return period corre- 
sponding to the double of a given flood is asymptotically the square of the origi- 
nal return period multiplied by e%. 


3. Estimation of the parameters for « = 0 


Since the graphical procedure is based on the logarithms of the droughts, it 
is appropriate to use them also for estimating the parameters. This means 
that, instead of estimating u, we estimate log u and 1/a' from the linear re- 
lation (2.7). The method of least squares can be used to minimize either the 
horizontal or vertical distances between theory and observations. The follow- 
ing compromise between the two methods, worked out for climatic extremes 
fil and for the linear relation between the extreme x and the reduced extreme 
y permits a rapid calculation and has been shown to be quite reliable. Accord- 
ing to this method, the parameter 1/a' is estimated by 


(3.1) s(log x) 


TN 


i. 
a' 


and log u is estimated from 
(3.2) log u = log x + Jy/a' 


This procedure requires the calculation of the mean log x and of the standard 
deviation s(log x) of the logarithms. The reduced mean Yn and the reduced 
standard deviation 0, depending only upon the sample size N are given in 
Table II. 
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TABLE 0 


ted means ¥y and standard deviations 7 of reduced extremes 


Expec 
N In oN N In oN 
15 5128 1.0206 48 5477 1.1574 
16 57 316 49 81 590 
17 81 411 -50 85 607 
18 5202 493 51 89 623 
19 20 565 52 93 638 
20 36 628 53 97 653 
21 52 696 54 5501 667 
22 68 154 55 04 681 
23 83 811 56 08 696 ‘ 
24 96 864 57 11 708 
25 5309 915 58 15 721 
26 20 961 59 18 734 
27 32 1.1004 60 21 147 
28 43 047 62 27 770 
29 53 086 64 33 793 
30 62 124 66 38 814 
31 11 159 68 43 834 
32 80 193 70 48 854 
33 88 226 72 52 873 
34 96 255 14 57 890 
35 5403 285 16 61 906 
36 10 313 78 65 923 
37 18 339 80 69 938 
38 24 363 82 72 953 
39 30 388 84 16 967 
40 36 413 86 80 980 
41 42 436 88 83 994 
42 48 458 90 86 1.2007 
43 53 480 92 89 020 
44 58 499 94 92 032 
45 63 519 96 95 044 
46 68 538 98 98 055 
47 73 557 100 .5600 065 
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Finally the theoretical values of the droughts are obtained from (2.7) as 


(3.3) log x = log u - y/a' 


If the probability paper is scaled in logarithms, the theoretical droughts can 
be immediately read off the straight line (3.3). In particular, we obtain the 
100 years drought by introducing the corresponding value of the reduced 
variate y = +4.6, and the expected drcught for any other desired return period 
by extending the line (3.3). 

It is interesting to compare two rivers such that 


> 


From (3.3) it follows for y = +4.6 that the first river which has larger dis- 
charges may have a more severe hundred years drought, i.e., a smaller value 
than the second river which has smaller discharges. This is not surprising 
since a river with a smaller modal flood may also have a higher 100 years 
flood than a second river with a higher modal flood, provided that the disper- 
sion for the first river is larger than for the second one. 

Table III shows the calculation of the parameters u and 1/a in the theo- 
retical line (3.3) for the Colorado River at Lees Ferry (Ariz.), 1922-39, 
(W.S.P. 879); Arkansas River, at Syracuse, Kansas, 1922-1939 (W.S.P. 877); 
Sta. Clara River near Central, Utah, 1911-30 (W.S.P. 879); Bourbeuse River 
at Union, Montana, 1922-39 (W.S.P. 877); Swift River at West Ware, Mass., 
1913-1945 (W.S.P. 1105); Colorado River at Glenwood Springs, 1904-09, 
1912-39 (W.S.P. 879); Chippewa River at Jini Falls Dam, Wisconsin, 1925-48 
and the Colorado River, Bright Angel Creek, 1924-39 (W.S.P. 879). The ob- 
servations and the results of the linear theory (3.3) are shown in graphs 1 
and 2. 


4. The general case 


The observed droughts traced in Graphs 1 and 2 are closely scattered 
about the theoretical lines (3.3). However, in other cases - see Graphs 3 and 
4 - the plotted points approach a curve which falls off less rapidly for in- 
creasing return periods and the assumption € = 0 is too strong. To account 
for these cases, we use the asymptotic probability P(x) that the smallest 
value exceeds x where the initial variate is limited by x 2¢ >0. From 
Table I 


(4.1) P(x) = exp [f+ | 


with 
(4.1') P(u) = 1/e; Pf) =1 


Thus u and € are location parameters and 1/a a scale parameter. The re- 
lation between mode and median is the same as in the special case. Therefore 
a pseudo-symmetrical case exists here too. 

For the graphical representation, the transformation (2.7) is used in the 
form 


(4.2) log (x -€) = log(u-€) -y/a' 


The observed droughts are again plotted on a decreasing scale on logarithmic 
extremal probability paper. However, the relation between the reduced variate 
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y and log x is no longer linear. The second derivative feet is negative. 
Therefore the curves log x, y are bent downward. As long as x is large com- 
pared to ¢, the curves are nearly linear. For increasing return periods T 
they approach a line parallel to the y axis at the asymptotic distance x = €. 
All curves with the same characteristic drought u intersect at the point x = u, 
y =0. For fixed values of u and ¢€ , the convergence toward the asymptotic 
value is more rapid as 1/a increases. Two curves log x, y are parallel if 
they have the same values of 1/a and u-€ although the parameters u and € 
themselves may differ. 

Since the previous graphical estimate of the parameters is no longer pos- 
sible, the classical method of moments is used. In the case a = 1, the proba- 
bility function (4.1) degenerates into an exponential function. The character- 
istic drought u becomes equal to the mean, and the lower limit is estimated 
from the observed mean X and standard deviation s of the droughts, as 


=X-s 


Here and in the following an estimated parameter is designated by the sym- 
bol *. 
In general we need 3 moments. The distribution p(x) of droughts obtained 


from (4.1) is 
a x-€\% 
px) = | 


Consequently the reduced moments (x- )K/(u- )K of order k are the gamma 
functions 


(4.3) (===) = 


u-€ 


For k = 1,2,3, three equations are obtained for the three parameters 
(44) r(1+2/a); (2=£) - = 7(1+34) 
The variance o? and the third central moment u, are 


(4.4') 
= (u- +2/a) +1/a)); = - - (x-€) + 


From (4.4) the skewness YB, defined by 


(4.5) 
becomes 
(4.5') 
VB, = (ra +8/a) - 3711 + + 1/a) + + 1/a)) (ra +2/a) 2 


This expression depends only upon the parameter 1/a. If the population value 
¥8, is replaced by the sample value Vb, , an estimate of 1/a is obtained. To 
facilitate this procedure, the right side of equation (4.5') as function of 1/a is 
given in Table IV, cols. 4 and 1.* The value of 1/a' in equation (4.2) is ob- 
tained from (2.7'). 


*This Table was calculated by Gladys R. Garabedian of Stanford University. 
I profit on this occasion to thank her for this important contribution. 
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The two remaining location parameters, the characteristic drought u and the 
lower limit «, are simple to estimate. Equation (4.4) for the mean leads with 
the help of (4.4') to 


(4.6) u=X+0 A(a) 


where 

(4.6) aa) = (1 - r(1+1/a)) (ria +2/a) - 

is given in Table IV, col. 3. Since 1/o has already been estimated from (4.5'), 
the parameter u may be estimated from (4.6) after replacing the population 
mean and standard deviation by the sample values. The result can be checked 
from the graph of the droughts on extremal probability paper with the help of 


the first equation (4.1'). 
To estimate the lower limit € , equation (4.4) for the first moment is 


written in the form 
_ X-ul(1+1/a) 
(4.7) (1 +1i/a) 
and the value of X is introduced from (4.6) and (4.6'). This leads, after trivial 
simplifications, to 


(4.8) €=u-¢0 B(@) 
where 
(4.8) B(a) = (r(1+2/a) - r1+1/a))-1/2 


is given in Table IV, col. 2. For estimating the lower limit € we use the pre- 
vious estimates of 1/a and u, and replace the population value ¢ in (4.8) by 
the sample value s. It is worthwhile mentioning that the lower limit is esti- 


mated here directly without successive approximations. 
Equations (4.6) and (4.8) lead to a criterion which decides whether or not 


zero should be chosen as the lower limit since 
(4.9) «20 if X+s (ACa) B(a) ) 20 
This condition may be written from (4.6') and (4.8') 
Ws 2 (1+1/a) - r(1+1/a)) 


whence 

(4.9") € 2 0 if%/x 

The lower limit is taken to be zero if the equality is fulfilled within the errors 
of random sampling. If € is negative and small, it can safely be assumed to be 
zero. The theory fails if the lower limit turns out to be larger than the ob- 
served smallest drought - provided this observation is reliable - and if € 


takes on a large negative value. 
In the pseudo-symmetrical case, 1/a = 0.30685 and the other two para- 


meters are estimated using (4.6) and (4.8) as 

fi = X + 0.34230s : € =u - 3.30649s 
The estimation of the three parameters is shown in detail in Table V for the 
Connecticut River, and in Table VI for other rivers. The three parameters 


are introduced into equation (4.2) which leads to log (x - €) and hence to the 
theoretical values log x and x as functions of the reduced variate y, i.e., of 


the return period T. 
439-10 


| 


TABLE IV 


Estimation of the Three Parameters 


Scale Multiple of Standard Deviations Reduced 
parameter for lower limit € for parameter u third moment 


B(a) A(a) VB, 
Equ. (4.8') Equ. (4.6') Equ. (4.5') 


78.981714 -448154 -1,081272 
39.989044 -446110 -1.024853 
26.986212 -443925 - .970702 
20.480808 -441603 - .918459 
16.574350 -439150 - .867967 


13.967343 -436568 -819101 
12.102862 -433863 -771740 
10.702446 -431038 -725772 
9.611395 -428096 -681102 
8.736889 -425043 -637637 


8.019861 -421881 -595296 
7.420934 -418614 -554002 
6.912848 -415245 -513687 
6.476131 -411778 -474287 
6.096505 -408216 -435743 


5.763261 -404563 -398002 
5.468210 -400822 -361012 
5.204984 -396996 -324729 
4.968556 -393087 -289108 
4.754903 -389100 -254110 


4.560770 -385036 -219697 
4.383495 -380900 -185835 
4.220878 -376693 -152490 
4.071085 -372419 -119634 
3.932577 -368079 -087237 


3.804052 -363678 -055272 
3.684400 -359218 - .023715 
3.572672 -354700 007458 
3.468048 -350129 038270 
3.369818 -345505 -068742 


3.277364 -340832 -098893 
3.190146 -336112 -128743 
3.107688 -331348 -158308 
3.029573 -326541 -187606 
2.955428 -321694 -216653 


2.884924 -316809 -245465 
2.817768 -311889 274055 
2.753697 -306935 -302437 
2.692475 -301949 -330625 
2.633890 -296935 -358632 


2.577752 -291893 -386468 
2.523887 -286825 -414147 
2.472138 -281734 -441678 
2.422364 -276622 -469072 
2.374435 -271490 496340 


2.328232 -266340 523491 
2.283647 -261174 -550535 
2.240583 -255993 -577481 
2.198946 -250801 -604336 
2.158655 -245597 -631111 


2.119632 -240384 -657812 
2.081807 -235163 -684448 
2.045114 -229937 -711026 
2.009492 -224706 -737553 
1.974885 -219472 -764038 


01 
.02 
.03 
04 
.05 
j 06 
.O7 
.08 
.09 
.10 
.12 
.14 
15 
-16 
-18 
19 
.20 
.22 
.23 
.24 
.25 
-26 
.27 
.28 
.29 
.30 
31 
.32 
.33 
.34 
.35 
.36 
37 
.38 
.39 
Al 
42 
43 
44 
45 
46 
47 
48 
.49 
.50 
51 
.53 
55 
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Scale , Multiple of Standard Deviations Reduced 
parameter for lower limit € for parameter u third moment 


B( a) A(a) 


1.941242 -214237 -790486 
1.908514 -209002 -816904 
1.876656 -203768 -843299 
1.845626 -19§537 -869677 
1.815385 -193311 -896045 


1.785897 -188090 -922408 
1.757128 -182875 -948772 
1.729045 -177669 -975143 
1.701620 -172473 1.001527 
1.674824 -167287 1.027929 


1.348631 -162113 1.054354 
1.623017 -156951 1.080808 
1.597958 -151804 1.107296 
1.573432 -146672 1.133822 
1.549420 -141557 1.160393 


1.525901 -136459 1.187012 
1.502857 -131379 1.213685 
1.480272 -126318 1.240415 
1.458130 -121278 1.267209 
1.436413 -116260 1.294070 


1.415109 -111263 1.321004 
1.394204 -106290 1.348013 
1.373683 -101340 1.375104 
1.353536 -096416 1.402280 
1.333750 -091517 1.429545 


1.314314 -086645 1.456904 
1.295217 -081799 1.484362 
1.276450 -076982 1.511921 
1.258002 -072194 1.539587 
1.239865 -067435 1.567363 


1.222031 -062706 1.595254 
1.204489 -058008 1.623264 
1.187234 -053341 1.651396 
1.170256 -048707 1.679655 
1.153550 -044105 1.708045 


1.137107 -039536 1.736570 
1.120922 -035002 1.765232 
1.104988 -030501 1.794038 
1.089299 -026035 1.822990 
1.073849 -021605 1.852093 


1.058632 -017211 1.881350 
1.043644 -012853 

1.028880 

1.014333 

1.0 


1. 2. 
-867491 2.309348 
-152233 2.640035 
-651524 2.996146 
-563330 3.382013 


-486053 -160077 3.802311 
-418382 -179747 4.262142 
-359209 -195656 4.767125 
-307573 -208071 5.323478 
-262625 217284 5.938118 


-223607 -223607 6.618761 
-038236 -191180 19.584859 
-005016 -115370 60.091733 
-000526 -062593 190.113240 


. 
1/a 
56 
57 
.58 
59 
.60 
61 
62 i 
63 
64 
65 
67 
.69 
1 
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TABLE V 


Estimate of the Three Parameters, 
Connecticut River at Sunderland, (Mass.) 


W.S.P. 1105 


Mean drought x from observation 1353.9 
Mean square x observation 1943682.143 
Variance s? from (1) and (2) 1105595.673 
Standard deviation from (3) 332.505 
Third power from (3) and (4) 36761607.56 
Mean cubic observation 2945659892. 


Third Central (1),(2),(6) 14673590. 
Moment 


Skewness Vb, = ms“ formula (4.5); 0.39915 
(7),(5) 


Parameter 1/a (8) Table IV, col.1 0.41458 
Parameter 1/@ formula (2.7') 0.18005 
Factor A(a) (9) Table IV, col.3 0.28957 


Char. drought a (1)(4)(11) formula 1450.2 
(4.6) 


Factor B(a) (9) Table IV, col. 2 2.55308 
Product sB(a)=@-@€ (4), (13) 848.9 
Lower Limit € (12),(14) 601.3 
Logarithm log (a -@) (12),(15) 2.92886 


Minimum Xap observation 780. 


18 Hundred years pn formula (4.2); 724. 
drought y=+4.6 
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TABLE VI 


The 3 Parameters for Different Rivers 


1 River Rhein Inn Arkansas R. Housatonic R. 
2 Station Rheinfelden Martinsbruck Salida (Colo.) Falls Village (Conn.) 
3 Period 1918-33; 35-1949 1917 - 1949 1912-18; '20-'39 1913-45 
4 Publication Schweiz. Hydrogr. Jahrbucher W.S.P. 877 W.S.P. 1105 
5 Number of years N 31 33 27 33 
6 Unit cubic meter per second c.f.s. c.f.s. 
7 Mean X 437.19 12.709 185.70 87.939 
8 Standard deviation s 77.620 1.50240 19.909 54.684 
9 Skewness yb, -60618 -27294 -24634 1.3660 
10 Parameter 1/@' -21309 -16052 -15636 -33730 
11 Char. drought @ 456.63 13.178 192.01 93.572 
12 Lower limit é€ 287.98 8.9406 134.58 18.078 
13 Observed Minimum xy 306.00 10.200 155.00 24.000 
14 100 years drought x,,. 302.7 9.484 145.90 20.000 


After calculation of the three parameters (lines 10-12) the theoretical 
droughts x obtained from (4.2) are plotted against y on logarithmic extremal 
probability paper (see graphs 3 and 4). The fit of the theory to the observa- 
tions is highly satisfactory, and the expected droughts for the desired return 
periods can easily be read off the graphs. 


5. Summary of procedures 


The extremal probability paper which is used by the Geological Survey 
[2,9] for floods can also be used for droughts if the logarithms of the droughts 
are plotted in decreasing order of magnitude. If the series of points (log x, y) 
is scattered about a straight line the limiting value of the droughts may be 
assumed to be zero. In this case, the two parameters are estimated from 
(3.1) and (3.2) and the theoretical droughts are obtained from (3.3). A criterion 
for the validity of the assumption ¢ = 0 is given in (4.9). 

If the series of points (log x, y) falls off less rapidly with increasing re- 
turn periods, the first three moments have to be calculated. Then the three 
parameters are estimated from equation:s (4.5') (4.6) (4.8) with the help of 
Table IV. The theoretical values of x calculated from equation (4.2) are 
traced as functions of y. 

If the number of observations is large enough so that the sampling errors 
are small and if it is reasonable to assume that the basic conditions prevail, 
the theoretical curves thus obtained can be used to estimate the most severe 
drought expected within a given number of years. Further, it is believed that 
this procedure may be useful for solving storage and irrigation problems. 
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PROCEEDINGS-SEPARATES 


The technical papers published in the past twelve months are presented below. Technical-division sponsorship is indicated 
by an abbreviation at the end of each Separate Number, the symbols referring to: Air Transport (AT), City Planning (CP), 
Construction (CO), Engineering Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
Sanitary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Water- 
ways (WW) divisions. For titles and order coupons, refer to the appropriate issue of “Civil Engineering” or write for a 
cumulative price list. 


VOLUME 79 (1953) 


JUNE: 194(CP) & (AT), 195(SM), 196(CP) & (AT), 197(HY), 198(ST), 199(EM), D-134(HY), D-141(HY). 


JULY: 200(SM)», 201(ST)», 202(EM)>, 203(sM)>, 204(AT)>, 205(EM)>, 206(sT)>, 207(8A)?, 208(8A)?, 209(ST)», 210(su)>, 
211(EM), 212(SU)», 214(HW)>, 215(SM)>, 216(ST)», 217(ST)>, 218(ST)P, 219(ST)>, 220(SM)>, 221 (HW)>, 
223(EM)”, 224(EM)?, 225(EM)?, 226(CO)”, 227(SM)?, 228(SM)”, 229(IR)?. 


AUGUST: 230(HY), 231(SA), 232(SA), 233(AT), 234(HW), 235(HW), 237(AT), 238(WW), 239(SA), 240(IR), 241(AT), 242(IR), 
243(ST), 244(ST), 245(ST), 246(ST), 247(SA), 248(SA), 249(ST), 250(EM)°, 251(ST), 252(SA), 253(AT), 254(HY), 255(AT), 
256(ST), 257(SA), 258(EM), 259(WW). 


SEPTEMBER: 260(AT), 261(EM), 262(SM), 263(ST), 264(WW), 265(ST), 266(ST), 267(SA), 268(CO), 269(CO), 270(CO), 271(SU), 
272(SA), 273(PO), 274(HY), 275(WW), 276(HW), 277(SU), 273(SU), 279(SA), 280(IR), 281(EM), 282(SU), 283(SA), 284(SU), 
285(CP), 286(EM), 287(EM), 288(SA), 289(CO). 


OCTOBER: 290(all Divs), 291(ST)°, 292(EM)°, 293(ST)°, 294(PO)°, 295(HY)°, 296(EM)°, 297(HY)°, 298(ST)°, 299(EM)°, 
300(EM)°, 301(SA)°, 302(SA)°, 303(SA)©, 304(CO)©, 305(SU)©, 306(ST)©, 307(SA)°, 308(PO)°, 309(SA)°, 310(SA)°, 311(SM)°, 
312(SA)©, 313(ST)©, 314(SA)©, 315(SM)°, 316(AT), 317(AT), 318(WW), 319(IR), 320(HW). 


NOVEMBER: 321(ST), 322(ST), 323(SM), 324(SM), 325(SM), 326(SM), 327(SM), 328(SM), 329(HW), 330(EM)°, 331(EM)°, 
332(EM)°, 333(EM)®, 334(EM), 335(SA), 336(SA), 337(SA), 338(SA), 339(SA), 340(SA), 341(SA), 342(CO), 343(ST), 344(ST), 
3$45(ST), 346(IR), 347(IR), 348(CO), 349(SM), 350(HW), 351(HW), 352(SA), 353(SU), 354(HY), 355(PO), 356(CO), 357(HW), 
350(HY). 


DECEMBER: 359(AT), 360(SM), 361(HY), 362(HY), 363(SM), 364(HY), 365(HY), 366(HY), 367(SU)®, 368(Ww)®, 369(IR), 
370(AT)®, 371(SM)®, 372(CO)®, 373(ST)®, 374(EM)®, 375(EM), 376(EM), 377(SA)®, 378(PO)®. 


VOLUME 80 (1954) 


JANUARY: 379(SM)®, 380(HY), 381(HY), 382(HY), 383(HY), 384(HY)®, 385(SM), 386(SM), 387(EM), 388(SA), 389(SU)®, 390(HY), 
391(IR)®, 392(SA), 393(SU), 394(AT), 395(SA)®, 396(EM)®, 397(ST)®. 


FEBRUARY: 398(IR)f, 399(SA)!, 400(Co)!, 401(8M)®, 402(aT)!, 403(aT)!, 404(mR)f, 405(PO)!, 406(AT)!, 407(sU)!, 
409(ww)f, 410(aT)!, 411(sa)!, 412(Po)!, 


MARCH: 414(ww)f, 415(su)f, 416(sm)‘, 417(sM)f, 418(aT)f, 419(Sa)f, 420(sa)f, 421(aT)f, 422(sa)!, 423(CP)f, 424(aT)f, 
425(sm)!, 426(IR)!, 


APRIL: 428(HY)®, 429(EM)®, 430(ST), 431(HY), 432(HY), 433(HY), 434(ST). 
MAY: 435(SM), 436(CP)°, 437(HY)®, 438(HY), 439(HY), 440(ST), 441(ST), 442(SA), 443(SA). 


a. Beginning with “Proceedings-Separate No. 200,” published in July, 1953, the papers were printed by the photo-offset method. 

b. Presented at the Miami Beach (Fla.) Convention of the Society in June, 1953. 

c. Presented at the New York (N.Y.) Convention of the Society in October, 1953. 

d. Beginning with “Proceedings-Separate No. 290,” published in October, 1953, an automatic distribution of papers was in- 
augurated, as outlined in “Civil Engineering,” June, 1953, page 66. 

e. Discussion of several papers, grouped by divisions. 

f. Presented at the Atlanta (Ga.) Convention of the Society in February, 1954. 
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